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Abstract 

We study Rayleigh-Benard convection in the high-Rayleigh-number and 
high-Prandtl- number regime, i.e., we consider a fluid in a container that is 
exposed to strong heating of the bottom and cooling of the top plate in the 
absence of inertia effects. While the dynamics in the bulk are characterized 
by a chaotic convective heat flow, the boundary layers at the horizontal con- 
tainer plates are essentially conducting and thus the fluid is motionless. Con- 
sequently, the average temperature exhibits a linear profile in the boundary 
layers. 

In this article, we rigorously investigate the average temperature and os- 
cillations in the boundary layer via local bounds on the temperature field. 
Moreover, we deduce that the temperature profile is indeed essentially linear 
close to the horizontal container plates. Our results are uniform in the system 
parameters (e.g. the Rayleigh number) up to logarithmic correction terms. An 
important tool in our analysis is a new Hardy-type estimate for the convecting 
velocity field, which can be used to control the fluid motion in the layer. The 
bounds on the temperature field are derived with the help of local maximal 
regularity estimates for convection-diffusion equations. 

1 Introduction 



1.1 Motivation 

Rayleigh-Benard convection is the flow of an incompressible Newtonian fluid 
in a container with an imposed temperature gradient due to heating of the 
bottom and cooling of the top plate. It is one of the classical models of fluid 
dynamics. With its numerous applications in geophysics, oceanography, me- 
teorology, astrophysics and engineering, Rayleigh-Benard convection played 
a central role in experimental and theoretical physics since the turn of the 
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last century. Depending on the heating rate, experimentahsts and numerical 
analysts observe a wide range of flow pattern: from purely conducting (i.e., 
motionless) states, over steady and oscillatory fluid motions, to chaotic pat- 
tern and fully developed turbulence. In view of the richness of the observed 
phenomena, it is not surprising that over many years, Rayleigh-Benard con- 
vection has become a paradigm for nonlinear dynamics, including instabilities, 
bifurcations and pattern formation, e.g. [U El [TT]. We refer to [TSl [131 El [I] 
for reviews and further references. 

We consider Rayleigh-Benard convection in the so-called high-Rayleigh- 
number and high-Prandtl-number regimes. This means, the applied temper- 
ature forcing is so strong that the fluid motion is unorganized. Moreover, 
inertia effects are negligible, so that, in particular, the observed flow pattern 
is rather chaotic than turbulent. Thinking of a container that is infinitely 
extended in horizontal directions, boundary effects on the vertical sidewalls 
of the container become irrelevant. In this situation, the flow pattern shows a 
clear spatial separation of the relevant heat transfer mechanisms: thin laminar 
boundary layers in the vicinity of the horizontal plates, in which heat is essen- 
tially conducted and a large bulk that is characterized by a convective heat 
flow. Since near the rigid walls of the container, the fluid is almost at rest, 
the relevant heat transport mechanism in the boundary layers is conduction, 
and thus, the temperature profile is essentially linear. Fluctuations around 
this profile take place on large length scales, with a small amplitude. It is in 
these boundary layers where the majority of the temperature drop between 
the hot bottom and the cold top boundaries happens. In the bulk, far from 
the rigid container walls, convection is only limited by viscous friction. The 
temperature is essentially equilibrated around the mean temperature of the 
system. However, due to the strong temperature forcing, overload heat is pro- 
duced on the bottom, which generates instabilities of the boundary layer: The 
boundary layer bursts and hot fluid parcels with warm tails, so-called plumes, 
detach and flow into the much colder bulk due to buoyancy forces. While 
rising, the plumes push aside the above fluid and are itself in turn deflected. 
In this way, the plumes take the forms of stalks with caps on their tops — 
the plumes obtain their characteristic mushroom-like shapes, cf. |20| 113]. (Of 
course, the corresponding effects can also be observed on the cold top plate. 
For convenience, we restrict our considerations to the bottom boundary.) 

The rigorous understanding of the dynamic flow pattern in the chaotic 
regime is far from being satisfactory. Until today, most of the work is devoted 
to the scaling of the Nusselt number, that is, the scaling of the average upward 
heat flux as a function of the imposed temperature forcing [6l O [71 [HI [T9l 116]. 
In this paper, we focus on the boundary layers at the horizontal bottom and 
top plates of the container. So far, boundary layer theories in Rayleigh- 
Benard convection were primarily derived as intermediate steps towards the 
understanding of the Nusselt number scaling. Since due to the absence of 
fluid motion, the boundary layers offer the main resistance for the heat flow 
through the container, the Nusselt number is dominantly determined by the 
boundary layer, and thus, for scaling theories, the width of these boundary 
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layer is of particular interest. (In fact, the Nusselt number is inversely propor- 
tional to the width of the layer.) However, because of their extreme thinness, 
experimental and numerical investigations of the boundary layers are very dif- 
ficult, and still, no generally accepted boundary layer theory is available for 
Rayleigh-Benard convection. An overview on the current state of experimen- 
tal and theoretical research for more general (than the one considered in the 
present paper) Rayleigh-Benard experiments can be found in [TJ Sect. VI]. 

The present work is a first attempt to partially characterize the observed 
pattern in chaotic high-Rayleigh-number convection — at least in the phe- 
nomenologically most regular region. Our aim is a rigorous justification of 
the laminar profile and thus the dominant role of conduction in the boundary 
layer. We rigorously establish local bounds on the temperature field and its 
gradients in the boundary layers that are (up to logarithmic corrections) uni- 
form in the system parameters (Theorem [1]) . These bounds indicate that the 
temperature field close to the horizontal plates is indeed essentially laminar 
and fluctuations only happen on relatively large length scales, with a weak 
dependence on the heating rate. In other words, our analysis proves that 
heat is transported essentially via conduction. Moreover, we can deduce that 
the average temperature decays linearly and with slope equal to the Nusselt 
number (Theorem [2]), as it is expected in conducting boundary layers. 

The remainder of the paper is organized as follows: In Subsection 11.21 we 
introduce the mathematical model and the Nusselt number; in Subsection II .31 
we present our main results and discuss the method of this paper. Section [2] 
is devoted to the analysis of the velocity field. Finally, Section [3] contains the 
proofs of Theorem [1] and [2l 

1.2 Model and Nusselt number 

Despite the complexity of the observed phenomena, the mathematical model 
for Rayleigh-Benard convection is relatively simple. If density variations of 
the fluid are sufficiently small and the thermal diffusivity is negligible com- 
pared to kinematic viscosity, the problem can be modelled by the infinite- 
Prandtl-number limit of the Boussinesq equations: 



Here, T is temperature, u the fiuid velocity, p the hydrodynamic pressure, 
and e the upward unit vector. We suppose that the container has the simple 
form [0,A)'^~^ x [0,//], where we refer to the first d — 1 coordinates as the 
horizontal ones, and to the last coordinate as the vertical one. We complete 
the system with periodic boundary conditions in the horizontal directions; at 
the rigid top and bottom plates, we set 



dtT + U- VT - I^T = 0, 
V-u = 0, 

— -|- = Te, 





(3) 




and u = for z e {0, H]. 



(4) 
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Thus, we suppose uniform heating/cooling of the bottom/top plate, and no- 
shp boundary conditions for the fluid velocity. We follow the convention to 
write the space coordinate and the velocity field as x = {y, z) G W^-'^ X E, and 
u = {v, w) S ]R^~-^ X R, respectively. 

The system is nondimensionalized and admits a single control parameter: 
the container height H. The side length A of the period cell is chosen arbitrary 
and has no significance in the subsequent investigation. We distinguish two 
regimes: 1) the linear regime for small container heights, <C 1. Here, 
heat transfer between bottom and top plate is exclusively due to conduction 
and the temperature field stabilizes in a linear profile. 2) the chaotic regime 
for large container heights, H ^ 1. In this regime, the fiuid experiences 
a strong temperature forcing, which leads to the formation of chaotic flow 
pattern, as described in the previous subsection. Notice that the present 
nondimensionalization differs from the common one in the physics literature, 
where the control parameter that measures the applied temperature forcing 
is the Rayleigh number Ra. Both, H and Ra are related via H = Ra^^^, 
which explains the term "High-Rayleigh-number convection" for the chaotic 
regime. For the sake of completeness, we like to mention that between these 
two regimes, there is a critical container height -ff* ~ 1 for which the linear 
profile becomes unstable and a steady circular fiuid (convection rolls) fiow 
sets in. These convection rolls become itself unstable for larger heights and a 
cascade of bifurcations can be observed, until the chaotic behavior sets in for 
1. 

The efficiency of the heat transport is measured by the Nusselt number 
Nu, which is defined as the average upward heat flux, 

Nu := ^ j\wT-d,T)dz i 1 j\wT)dz + ^. 

Here, the brackets ( • ) denote the horizontal space and time average, i.e., 

1 1 



1 r 

if) ■= limsup— / 

to too 1^0 JO 



f{t,y)dydt 

[0,A)<^-i 



for any function f{t,y). Some equivalent expressions for the Nusselt number 
can be derived: 

Nu = {wT - d,T) for all z G [0, H] (5) 
|VT|2)dz (6) 







(7) 



Indeed, for ([5]), we average ([T]) in horizontal space and time and apply the 
boundary conditions (jl]). Identity ([6]) can be derived by testing the heat 
equation ([T|) with T—{l — z/H), integrating by parts and using the divergence- 
free condition ([2]) and the boundary conditions ([4|), while ([7]) follows from 
testing ([3]) with u, integrating by parts and using ([2]) and ([H) again. 
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Over many years, the dependence of the Nusselt number on the control 
parameter H has been studied extensively, cf. [121 [1] and references therein. 
In the absence of inertia, as in the present model, the expected scaling of Nu is 
the following: In the linear regime H <^ 1, when heat transport is essentially 
due to conduction, the Nusselt number scales like the imposed temperature 
gradient: Nu ^ 1/H. In the chaotic regime H ^ 1, it is conjectured that 
conductive and convective heat transport are balanced in the system, in the 
sense that Nu ~ 1. While the first scaling law can be easily established, cf. 
[n\ Theorem 2], a rigorous derivation of the second one remains reluctant 
until today. The presently best upper bound on the Nusselt number in the 
chaotic regime was derived by Otto and the author in [16] and is optimal only 
up to a double-logarithmic factor: 

Nu < ln^/'\lnH) as if > 1. (8) 

Observe that proving rigorous lower bounds on Nu is substantially different 
from proving upper bounds: Lower bounds depend sensitively on the partic- 
ular choice of the initial data. Indeed, there are ungeneric solutions to the 
system for which the heat flux is less efficient, for instance the purely 

conducting state T = 1 — z/H, u = 0. In this case, it is Nu = 1/H ^ 1. 
Therefore, one can only expect to prove (physically relevant) a-priori upper 
bounds on the Nusselt number in the chaotic regime. 

1.3 Main results of this paper 

From now on, we restrict our attention to the chaotic regime, that means, we 
assume that 

if » 1. 

Moreover, in order to rule out "ungeneric" configurations like purely conduct- 
ing solutions (as described in the last subsection) or steady fluid motions, e.g. 
convection rolls, we focus on solutions for which 

Nu > 1 as ii > 1. 

We are interested in statistical properties of the temperature field in the 
boundary layer in the chaotic regime. Our first result gives upper bounds on 
derivatives of the temperature field up to fourth order. We show that in a 
boundary layer of order one, the temperature can be controlled in appropriate 
Sobolev norms uniformly in terms of H — modulo logarithmic pref actors. 

Theorem 1. Assume that H :$> 1 and Nu > 1. Let a G N, q < 4. Then 
there exists a 7 > such that 

/\|V°r|^/")dz < In^'H. (9) 
Jo 

In the above statement, the size of the boundary layer is not explicitly 
fixed to one. Instead, we can rather think of these estimates to hold in any 
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boundary layer whose width is of order one or logarithmic in H. However, for 
simphcity we restrict ourselves in the following to the order-one case. Likewise, 
for symmetry reasons, the same bounds hold in the upper boundary layer. 
Moreover, for the benefit of a concise statement and to trim the estimates in 
the proofs, we do not compute explicit exponents 7 on the logarithms. Upper 
bounds on 7 (at least for the second and third order derivatives), can be found 
in the author's PhD thesis |17j . 

The bounds stated in Theorem [1] can be considered in the context of 
rigorous upper bounds on the Nusselt number initiated by Constantin and 
Doering. Indeed, in view of ([6]), the bounds in [5l[8l[16], e.g. ([8|), can be read 



for some 7 > 0, and thus, estimate ([9]) appears to be a new contribution in the 
study of universal bounds on the temperature field in infinite-Prandtl-number 
Rayleigh-Benard convection. 

The estimates in Theorem [1] can be used to show that, in the vertical 
boundary layers, the temperature profile is essentially linear: 

Theorem 2. Assume that H ^ 1 and Nu > 1. There exists a 7 > such 
that for any z € [0, 1] 



Again, by symmetry, an analogue statement can be derived in the upper 
boundary layer. As the majority of the temperature drop between the bottom 
and the top plate in the Rayleigh-Benard experiment just happens in the 
boundary layers, the heat flow must be inversely proportional to the width of 
the layer. The precise slope can be easily deduced from ^ by choosing z = 
and exploiting the boundary conditions (HD for w, namely dz\^^Q {T) = —Nu. 
Now, in Theorem [2l we recover this slope all over the boundary layer. The 
main benefit of the above estimate is the cubic control of the deviation in 
vertical direction around the linear profile, which only depends logarithmically 
and thus weakly on H. 

The method we use in this paper in order to derive the bounds on the 
temperature field stated in Theorem [1] are motivated by bounds derived by 
Otto for the viscous Burgers equation in [T3] and Otto and Ramos for the 
Navier-Stokes equation in |15j . In both papers, the authors establish Lp 
maximal-regularity-type estimates with tools borrowed from Harmonic Anal- 
ysis. Using and refining the techniques from these two papers, we derive 
new local maximal regularity estimates for convection-diffusion equations, see 
Propositions [2] and [3] on p. [TSF. After differentiating the heat equation ([1]) in 
spatial directions and and localizing in the boundary layer, these estimates 



as 




\{T\^)-{l-zNu)\ < zHrCH. 



It follows from the above Theorem that 




1 - zNu for z < \ir^/^ H. 
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can be applied to gain control on gradients of the temperature field. We have 
to comment on the choice of the norms in the assertions of Theorem [TJ These 
are determined by the leading order error terms, which come from differen- 
tiating and localizing the transport nonlinearity. Only few tools are known 
that are applicable in order to bound these terms uniformly in H (modulo 
logarithms): localized L2 bounds on the velocity field in the spirit of those 
derived in [H], cf. Proposition [T] on pagell2t the maximum principle for T, and 
interpolation inequalities of the form 

j |v"-ir|'^/(°-i) dx < sup |r|^/"("-i) j |v"r|^/" dx. 

The value of p in the Lp maximal-regularity estimates is chosen optimally 
with respect to the estimates at hand. 

2 Bounds on the velocity field 

In this section, we collect a series of L2 bounds on the velocity fields: some 
known results from ^ and [16] and some new results together with their 
proofs. The main new result, stated in Proposition [1] below, is an estimate 
on the shear velocity v in the boundary layer. As this estimate might be of 
independent interest, in this section, we state all results for boundary layers 
of arbitrary width 5. 

We consider the stationary Stokes equation for the velocity field u = 
{v, w) G W^^^ X R and the hydrodynamic pressure p in the Rayleigh-Benard 
problem, i.e., equations (H]) & (l3|) equipped with periodic boundary conditions 
in all d — 1 horizontal directions and with no-slip conditions on the vertical 
boundaries, i.e., w = v = for z € {0,H}. We observe that, because of ([2]), 
we have the additional information that dzW = — Vy ■ v = for z G {0,H}. 
It is convenient to eliminate the pressure term in ^ with the help of the 
incompressibility condition This leads to a fourth-order equation for w: 

A^w = -AyT and w = d^w = for z e {0, H}. (10) 

Here, Ay denotes the Laplace operator in the horizontal variables. Indeed, 
applying the divergence operator to Q yields Ap = dzT, and therefore, 
applying the Laplace operator to the vertical component of ([Sj), we obtain 
A'^w = AdzP — AT = —AyT. In a similar way, making additionally use of 
the boundary conditions, we see that v satisfies the equation AyV = ^y^y -v, 
and thus by ([2]), u is determined by 

V = {-AyY^VydzW. (11) 

Furthermore, since there are no external forces acting on the fluid velocity, 
we have 

(n(to)) = 

for every to > 0, which follows immediately from averaging ([2]) and ([3]) in 
horizontal direction, and exploiting the boundary conditions for u. 
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At this point, we like to recall some bounds on the vertical velocity compo- 
nent w, derived in [8] and [16]. The key estimate of [8] (here, slightly extended 
by including the horizontal gradient |Vyw|^) is the following: Let w and T be 
periodic in y and satisfy (jlOp . Then 

"^i,<fm,,. (12) 

In Lemma [3] below, we will extend this result to the horizontal velocity com- 
ponent V. Notice that, as a byproduct of the result in [8], the r.h.s. of (fT2|) is 
controlled by the Nusselt number (modulo a logarithmic correction): 

/^i^dz < {\nH)Nu. (13) 
Jo z 

We also refer to Lemma [161 Lemma 2] for a direct proof of this estimate. 
Consequently, ^ k ([HD imply 

Lemma 1 ([8l|16]). 

/ ^' ' ' dz < {In H)Nu. (14) 
Jo z 

In Lemma 3, see also eq. (18)], the authors derive an L2 maximal 
regularity estimate for the third-order derivatives of w from (|lUp . By the 
Nusselt number representations (l6|)&:([7]), the result can be stated as follows: 



Jo 



(15) 



or, invoking the formula (jlip : 
Lemma 2 (|16]). 

rH 

{\V^u\'^)dz < Nu. (16) 



The H2 ^ norm in (jlSp can best be understood on the Fourier level: For 
any periodic function f{y) we define 

Here J-f{k) is the Fourier transform of / at wave number A;, i.e., 

1 



^fik) := [ e^'-'f{y)dy. (17) 



Our first new result is the following H2 ^ bound on dzW, which extends 
estimate (|12p to the horizontal velocity component v via the relation (jlip : 
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Lemma 3. Let w and T be periodic in y and satisfy (|10|) . Then for any 
< 5 <^ H and a > 1 it holds 







and 



(18) 

" + r ^^^1^ < ff ) ... (19) 

^ J(5 ^ \ J Jo z 



As we shall see, the main ingredients for estimates (jlSp and ()19p are Hardy- 
type inequalities for the function i.e., inequalities of the form 

^1 d 

^ / ^'-^i-jfdz (20) 



for some G R. This inequality fails for = — but only logarithmically. 
We exploit this logarithmic failure in two different ways. In the "bulk" (6, H), 
the logarithmic failure of (j20p with u = produces the prefactor ln'^{H/6) in 
(jl9p . Notice that the estimates leading to (jl9p are sharp. In particular, the 
logarithmic prefactor can be recovered in the critical Hardy inequality when 
choosing (f> w ln{z/5) (with the appropriate boundary conditions at z = H). 
In the "boundary layer" (0,(5), we prove a subcritical Hardy inequality, i.e., 
(j20p with some u = —a < 0, which leads to (jl8p . Notice that in view of the 
prescribed boundary conditions for w, it holds ^l^a/^^."'! ^ ^ (|V,; ^8z-u)| ) ^ ^ 
for z <C 1, which just beats the integrability. Consequently, we can expect 
(j20p to be true only for u < 0, i.e., a > 0. In light of this observation, our 
bound for the boundary layer term seems to be suboptimal, since it requires 
a > 1. However, it is optimal in terms of scaling — a fact that will be 
exploited in Proposition [1] below. The logarithmic prefactor in (jlSp stems — 
roughly speaking — from the fact that the boundary terms at z = 5 have to 
be estimated via the bulk estimate (fT9]). 



Proof of Lemma O We prove the result on the Fourier level, cf . (|17p . Under 
Fourier transformation, (jlOp translates into 



\k\^J^w — 2—-^Tw + -—-pr—-^Tw = FT 
dz"^ dz^ 

and 

= —Tw = for z G {0, H}. 

dz 

We start by recalling the following inequality that has been derived in [HI 
p. 238&;239] and that can be easily reproduced by testing the equation for Tw 
with z~^T'w and integrating by parts: With <^ = ^ we have 

J_^",3|^^|.,, < ^"i^TjSd.. (21) 
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Notice that we may assume — applying an approximation argument — that 
(j) has the same boundary values as J-w, so that boundary terms vanish when 
integrating by parts. Now the statements in (|18p and ()19p follow from the 
three Hardy(-type) estimates 

"^Wdz < W(!L\ f z\^^4'?dz (23) 



and 



LFu; dz + I — — + / -^\-r-'w\ dz 



Indeed, for instance estimate (1181) can be derived as follows: 



\J-W\ dz + -r—TT / ^ \—J-W\ dz 



|A;|2 z5-al I |^|2 ^3-al^2 

9 1 /-^ , d .2 , 1 1 



< In^ ( ^ 1 / -TwTTdz. 



In view of the Plancherel Theorem, summing over all wave numbers k and 
averaging in time yields (|18|) . Estimate p9|) is derived similarly. 

It remains to prove statements (j22p - (j24p . Observe that (j22p is a standard 
Hardy estimate. We omit its proof which is very similar to the one of the 
critical Hardy inequality (j23p . The latter follows from 

< (In^i?) ^l^'Al'd^ 

via the transformation z = 6z, H = SH, and (j){z) = 4'{z). This estimate is 
easily established. We integrate by parts and estimate with the help of the 
Cauchy-Schwarz inequality: 

^ 1 ~ d 

-\^\^dz = / —{InzM^dz 
z Ji dz 

< 2 [""{Inzm^Jldz 



< 2{ln H) I j^" ^14)1"^ dz J" z\-^4>\'^dz 



1/2 
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which immediately yields 



The proof of (I24p is slightly more involved. We introduce the weight 
function 



and have the obvious estimate 



4^ for < 2 < 5, 
z for 5 < z < H, 



p{z)\—<j)\^ dz < / z|-^<^|2dz, (25) 







dz /n dz 



because of p{z) < z. Varying the weight is helpful, since it allows to treat 
both boundary layer and bulk estimates simultaneously, and thus avoids the 
appearance of boundary terms at z = (5 (for which we do not know an appro- 
priate estimate). Using the definition of (/>, we easily compute 

\ p{z)\—4)^ dz = j p{z)—x\—J^w\'^dz — 4Re f p{z)—^T'w—Fwdz 
Q dz Jq z dz Jq z dz 

+ 4/ p{z)^\Fw\'' dz, 
Jo ^ 

where Re Z denotes the real part of a complex number Z. We see via inte- 
gration by parts that 



S' 



— = — \J"w\ dz 



a 







dz \ z 
f-H 1 rH 



f 1 fid 

: / p(z)—^\J'wfdz — 2Re / p{z)—^Fw—Fwdz, 
Jo Jo ^ dz 



so that we may rewrite the above formula as 
p{z)\^cP\^ dz 



dz 

1 I . d 



—Tw\'^ dz 



2 ^ 1 ^ ~d~r , a /"^ 1 ,^ ,2 , 

— —Re / — I — J^w—j^wdz + — / — = — \J^w\ dz 

1 d 1 ~d 

+ / -^\—Twfdz-2Re / —rTw—Twdz. (26) 
J& z-^ dz Js z'^ dz 

We invoke the Cauchy-Schwarz and Young inequalities to deduce 
^" p{z)\^^(t>?dz 

a-1 1 , ^ ,2 , 1 /"^ 1 , ,2 , /"^ 1 , T- ,2 , 

> — ; — / — E — \J^w\ dz ^ — / —:t\—J^w\ dz — 2 —^IJ^wl dz, 
- 6^ Jo ^s-"! I 2Js z^Uz ' Js z^ 
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which becomes 



— 1^ — / -5 — \-rw\ dz + I -i;\—J-w\ az 
5^ Jo ' Js z^Uz ' 



in view of the definition of (j) and formula (j25p . Now, we obtain control on 
the second term on the l.h.s. of (I24p via ()26p and the Young inequality. This 
completes the proof of Lemma [3l □ 

The following proposition is a consequence of the Lemmas [H [2] & [3l 

Proposition 1. For any < 6 < H it holds 

/ {\Vyv\'^)dz+ / {\Vw\'^)dz < 6{\nH)Nu. 
Jo Jo 



Moreover, if jj <^ 6 ^ In H then 

f 

Jo 



{\d,v\'^)dz < 6{\n^H)Nu. 



As an immediate consequence of the previous proposition together with the 
Nusselt number bound ([8]), we observe for later citations that, in a boundary 
layer of order one, we have the estimates 

[\\u\')dz < f\\Vu\^)dz < In-yR, (27) 
Jo Jo 

for some 7 > 0. Here the first estimate is due to Poincare's inequality. 

Proof of PropositionUl The bounds on Vu) and VyV, follow directly from 
Lemma [T] via 



/ {\Ww\'^)dz < 6 / ^ ^dz < 5(\nH)Nu 

Jo Jo ^ 



{\Vyv\^)dz < / {\V-^VydM^)dz = / {{d,wf)dz. 



and 

/ {\yyV\')dz < / (|V, V, 

Jo Jo Jo 

It remains to estimate the dzV-terui. Because of (jlip . we have to show 

/ {\Vy'd^M^)dz < 6{ln^H)Nu. (28) 
Jo 

First, we apply the elementary estimate 

12 



to / = V y ^dzW and obtain 

S 5 S 

^ {\V~'d'M')dz < ^J^ {\Vy'dM')dz + S^ {\V-'dtw\')dz. 

For the first term on the r. h. s. of the above equation, we invoke the boundary 
layer estimate from Lemma [3] with a = 3, 

< dHn^ (j\ (In H)Nu 

< 5^{ln^H)Nu. 

The second term is estimated with the help of Lemma 3 in [16] : 

{\V-'dtw\^)dz < Nu. 
We combine the above estimates and deduce (1281) since 5 <C Inff. □ 



We conclude this section with a global bound on the velocity field, which 
combines classical Calderon-Zygmund theory with the maximum principle for 
the temperature. 

Lemma 4. For any 1 < q < oo, it holds 

±-l\\Vundz < j\\V^u\^)dz < H. (29) 

Proof of Lemma [7} The first estimate in ()29p follows immediately from the 
Poincare inequality and the boundary conditions The second statement 
follows from Lq maximal regularity estimates for the Stokes equation. 



''{\V\\'^)dz < ["{{rndz, (30) 
'0 Jo 

and the maximum principle for the temperature in the sense that 

(sup|r|«) < 1, 

X 

cf. ()4ip below and discussion on page [T71 Although ()30p is a very classical 
estimate, we have to argue that the estimate is uniform in the aspect ratio ^. 
By rescaling x, u, and p, and redefining A, we may w.l.o.g. assume that H = 1. 
Furthermore, since the pressure is unique up to a constant, we suppose that 

/ pdydz = 0. (31) 
i(o,2)''-i 
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Our starting point is the analogue result for the Stokes equation on an infinite 
strip. Consider 

-Au + Vp = /, (32) 
V-u = g (33) 

m M'^-i X (0, 1) with non-slip boundary conditions at the vertical boundary: 
u = for z € {0, 1}. Then we have 

/ \V^u\'' dydz + [ [ \Vp\'^dydz 
< [ [ \f\'^dydz+ [' [ \Vg\Uydz (34) 

Jo Jud-^ Jo J]R<*-i 

by classical theory, cf. \10\ Chapter IV], as long as the r.h.s. is finite. In order 
to apply this estimate, we introduce a smooth cut-off function t] = r]{y) with 
the properties that 77 € [0, 1], 77 = 1 for \y\oo < NA, = for |y|oo > NA + 1, 
and sup^ \^yV\ ^ 1) sup^^ l^yVl ^ 1- Here, € N is an arbitrary, but fixed, 
number that has to be chosen explicitly at the end of the proof. It is readily 
checked that 

u = r]u, p = 'i]p, f = r]Te- Ayriu-2Vyi]-VyU+ ^^^^p, g = Vyi]-v 

solve ([32|) &: (j33|) . In particular, if we use the cut-off properties of r], we deduce 
form (j34p the estimate 

[ [ \V\\Uydz+ [ [ \Vp\'^dydz 

Jo J\y\^<NA Jo J\y\oo<NK 

< /V \T\''dydz 
Jo J\y\oo<NA+l 

+ 11 {\u\'' + \Vu\'' + \p\'')dydz. 

Jo J N K<\y\^<N \+l 

In view of the periodicity of T, n, and p in y, we may rewrite the above 
equation as 

N'^-^ C [ dydz + N'^^^ C [ \\/p\Uydz 

Jo J[0,A)rf-i Jo J[0,A)rf-i 

< N'^-^ [ [ \T\idydz+ [ [ \T\idydz 

Jo J[0,A)'*-i Jo J(0,2)rf-i 

+ C I {\u\'i + \Vu\'^ + \p\'^)dydz. 

Jo J(0,2)d~i 

Notice that, by the homogeneous boundary conditions of u at z G {0, H}, we 
can apply the Poincare inequality in the vertical variable both for u and 
and obtain 

/ / \u\'^dydz < [ [ iVulUydz < [ [ \V\\Uydz. 

Jo J{0,2Y-^ Jo J{0,2Y-^ Jo J{0,2Y~'^ 
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Moreover, thanks to the Poincare inequahty in y yields 
/ \p\'idydz < C ! \Vp\'^dydz. 

Jo J{0,2)'i-^ Jo J(0,2)d-i 

We conclude that 

r-l r rl 



N"^-^ [ [ \V^u\Uydz + N"^-^ [ [ \Vp\Uydz 

Jo i[0,A)'*-i Jo i[0,A)'*-i 

rl r rl 

rd-l 



< N"-' / / \T\Uydz+ / / \T\Uydz 

Jo i[0,A)'i-i Jo i(0,2)d-i 

+ / / \V^u\'^dydz+ [ [ \Vp\'^dydz. 

Jo i(0,2)'*-i Jo J(0,2)<*-1 

Thus, choosing N sufficiently large (A^A > 1 if A ^ 1), we may absorb 
velocity and pressure terms into the 1. h. s. Normalizing by the system size 
and averaging in time, this proves ()30p . □ 



3 Estimates on the temperature field 

In this section, we prove Theorem [1] and [2l We first show, how Theorem [1] 
implies Theorem [2l 

Proof of Theorem\^ The proof is elementary. Observe first, that 

{d%=,T)\ < In^H. (35) 

Indeed, this is a direct consequence of the Poincare and Holder inequality and 
([9]) with a = 3 and a = 4: 

-1 rl 



{d%=oT) < [ {\dlT\)dz+f {\d%,=,T-dlT\)dz 
Jo Jo 

{\dlT\'/') dz^ + {\dtT\)dz 



< InT H. 



By Tailor expansion, we may write the average temperature as 

m = {T\,=o) + z{d,\.=oT) + \z'{d%^^T) 



+ \zHd%^,T) + ^J\dtT){z-~zfd~z 



We first identify the expressions on the right up to the quadratic term with the 
help of the boundary conditions ([4]): We immediately see that (Tj^^o) = 1) the 
Nusselt number representation ([5]) evaluated at z = yields {dz\z=oT) = —Nu, 
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and exploiting (HD in ([T]) again, we see that (3^|^^qT) = 0. Thus, the above 
identity turns into 

{T^,)-il-zNu) = U^{d%^^^T) + ^^j\dtT){z-~zfd~z. 

To bound the terms on the r.h.s., we invoke ([35]) and estimate 

rz j-l 1(9} 

/ {dtT){z - zf d~z < / {\dtT\)dz < z'lrPH. 
Jo Jo 

This proves Theorem [2l □ 

For the proof of Theorem [H we need some preparations. We start with 
three interpolation inequalities, stated in the Lemmas [5] and [6l 

Lemma 5. Let denote a periodic function with vertical period 2H . Then 

pH N 3/4 . H X 1/4 / X 1/2 

J jVC\'/')dzj <i^J^{e)dzj ^{\V'C\)dzj . (36) 



Proof of Lemma O We first observe that it is enough to prove estimate ([36]) 
in the case where ( = C{s) is a periodic function on some interval [0,a], i.e., 

f-a \ 3/4 / .a \ 1/4 / .a \ 1/2 

/IC'I-..) <(/c^..) (/,C"M») . 

The general case follows (roughly speaking) when choosing s = Xi, integrating 
over the remaining d — 1 variables, and summing over all i. Moreover, the 
above statement is a consequence of the estimate 

i-b \ 3/4 y \ 1/2 pb 

J \C'\^^''dsj < i^J C^dsj +J \C"\ds for ah 6 >0, (37) 

as one can see by rescaling length by a/b and optimizing in b. We now turn to 
the proof of ([57]) . If 5 <C 1, then ([57]) reduces to a Poincare-type inequality. 
Indeed, making use of the periodic boundary conditions of we estimate 

^ pb \ 3/4 ^ pb 

/ IC'l^/'ds < 6'/'supjC'| < fc^/M \C"\ds < / \C"\ds. 
kJo J s Jo Jo 

Observe that ([37p is a standard Ehrling estimate in the case where 5 ~ 1 — 
regardless of the actual boundary conditions. Therefore, if 6 ^ 1, we can 
invoke a summation argument based on the "6 ~ l"-case. For instance, if 
6 = IT G N, we have 

rK \ 3/4 K , \ 3/4 



V \Jk-l J Jk-1 . 



f j-K \ 1/2 r-K 

/ Q^ds] + \C"\ds. 
Jo / Jo 

This proves Lemma [5] □ 
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The following estimates are designed to make use of the maximum principle 
for the temperature field \T\ < 1, cf. ()4ip below. 

Lemma 6. Let C denote a periodic function with vertical period 2H . Then 
r (|VC|'> dz < sup ICP r (IV^CP) dz (38) 

and 

^" (I V^CP) dz < sup |Cp/=^ {\V'C\^/') dz. (39) 

Proof of Lemma\^ We start with estimate (j38p and fix i G {1, . . . , d}. Inte- 
gration by parts and applying the Holder inequality yields 

mf)dz = -3 r {adiCfd^Odz 

-H J~H 



< 



sup Id ( / ^mo')dz^ i^j ^{{dHf)dz 



It remains to apply the Young inequality and sum over all i to deduce (j38p . 

Now, we turn to ([3SD- Again, it is enough to show this statement for 
V^C replaced by dfC- The estimate for mixed second derivatives follows via 
integration by parts and the Cauchy-Schwarz inequality. The proof is very 
similar to the one of (j38|) . We integrate by parts and use the Holder inequality: 

{idfC?)dz = - r {d,(:dK)dz 
-H J-H 

< i^J ^m)')dzj i^J JdfC\'/'')dz 

In order to obtain (j39|) . we just insert (p8|) (with V^C replaced by df( as in 
the derivation of ([55]) ) in the above estimate. □ 

To study the convection-diffusion equation ([T]), it is convenient to rewrite 
this equation in terms of ^ = T — (1 — z/H), that is, 

dt9 + u-Ve - AO = ^w. (40) 
H 

This is beneficial since has homogeneous boundary conditions: ^ = for 
z G {0,H}. We may w.l.o.g. assume that 

sup|r| = sup|6'| < 1. (41) 

X X 

In fact, by the maximum principle for the temperature, T € [0, 1] (and thus 
9 € [—1,1]) is preserved during the evolution. Even if this condition is not 
satisfied initially, T attains this temperature range exponentially fast. Since 
we are working with long time averages, finite time intervals in which (j4ip is 
violated will not affect our results. 



17 



Notice that by the definition of 9 and the Nusselt number representation 
I, we have 

rH 

{\Ve\^)dz ~ Nu. (42) 

For higher order derivatives, it holds V^T = V^0, V^T = and so on. 

We observe that odd reflection of 9 and w and even reflection of v leaves 
the equations ([KJ]) and ([2]) invariant. Regarding the boundary conditions of 
0, w, and V, we infer from ()40p that = for z S {0, i/}. Hence we may 
think from here on of as a (restriction of a) 2i7-periodic function in 
z, of w; as a (restriction of a) 2//-periodic function in z, and of f as a 
(restriction of a) 2i7-periodic vector field in z. In particular, 9 is four 
times weakly differentiable, w is twice weakly differentiable, and v is one time 
weakly differentiable in z. 

The core of the proof of Theorem [1] are the following maximal regularity- 
type estimates for convection-diffusion equations: 

Proposition 2. Let C,, u, and f be smooth periodic functions with vertical 
period 2H satisfying 

dtC + u-VC-AC = /, (43) 
V-n = 0. (44) 

Assume that u = for z = 0. Assume further that 

supple C [-5,5] (45) 

for some < 6 < H . Let 1 < r < 2. Then 

H \ 
{\V^Q\'')dz 

-H 

< Mi^j jfndz 

H \ lA 



+ e-^'l^j ^{\VfY)dz^ 



H X 1/2 

(C^) dz 
-H J 

25 \ 1/2 . .H \ {2-r)/2r 

{\vu\')dzj jcr^^'-^-^)dzj 

H \^/'^ / f-H \ (4-r)/4r 



+iv(^y" ^{\vu\^)dz^ (y"^(icp"/^'""^)d^' 

+ e-^ n ^(|Vtx|^) dz\ U ^(|VC|"''/('-^-)) dzj (46) 



for any M,N G Nq. 

In the case r = 2 we can do better: 
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Proposition 3. Let C,, u, and f be smooth periodic functions with vertical 
period 2H satisfying 

dtC + u-VC-AC = f, (47) 
V-u = 0. (48) 

Assume that 

supp.Cc [-,5,5] (49) 

for some < 6 < H . Then 

r (I V^CP) dz < r dz + sup Id {\Vu\^) dz. (50) 

J-H J-H X J-S 

In the proof of Theorem [H we will apply Propositions [2] and [3] to the local- 
ized temperature field and its derivatives, i.e., C = di{i]6), dij{i]6), where 
r] = i](z) denotes a cut-off function. Before, we like to comment on ()46p and 
(j50p . The derivation of a maximal regularity estimate for convection-diffusion 
equations in L2 is a straight-forward computation, cf. proof of Proposition 
[3l Its Lr analogue, however, is challenging. The difficulty in deriving scale- 
invariant Lr estimates stems from the fact that, for our purpose, the con- 
vection term u ■ VC can not be treated as a error term: In fact, local (and 
thus, up to logarithms uniform) bounds on the velocity field are known only 
in L2, see the previous section. Hence ordinary maximal regularity theory in 
Sobolev spaces is not suitable for our purposes. On the level of Besov norms, 
however, the transport-nonlinearity can be controlled in an elegant way. At 
first glance, we expect an estimate of the type 




Compared to this estimate, the additional error terms that occur in (|46p are 
caused the fact that Sobolev and Besov norms are in general not equivalent. 
These error terms produce an additional logarithmic prefactor in Theorem [TJ 

We postpone the proof of Propositions [2] and [3] to the end of this paper 
and start with a series of global bounds on the temperature field, all based on 
Proposition [T] 
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Lemma 7. There exists a 7 > such that 



r{\VT\')dz < ["{\V'T\^)dz 
Jo Jo 

< r{\V'Tf^)dz 
Jo 



< 



H 

4 



(|v^r|)dz 







< {In-r H)H. (51) 

Proof of Lemma We first remark that it is enough to prove the statement 
with T replaced by 9. Observe that in view of the Lemmas [5] and [6] above, 
which we apply to ^ = and C, = 9, respectively, and the maximum 

principle for the temperature ([^T|) . it is 

r {\V9f)dz < (" {\VH\^)dz, (52) 

J~H J~H 

{\VH\^)dz < r {\vH\'"'')dz, (53) 

-H J-H 

r {\VH\'l^)dz < r {\V^9\)dz. (54) 

Indeed, ([52]) and ([53]) follow immediately from ([38]) and ([39]) . respectively, 
with C = 6* via (dH). For ([MD, we apply ([36|) with C = V^6' and estimate: 



H 13611 / /•« \ ' / 

^{\VH\"'')dz < i^J ^{\V'9\^)dzj i^J jV^9\)dz 

^ / rH \ 1/3 / N 2/3 

< i^J jV'9\^/')dzj i^J jV'9\)dzj 



It remains to invoke Young's inequality to obtain ([54 

As a consequence of ([52]) - ([M]) . for ([5T]) it is enough to prove 



^0 



dz < H)H. (55) 



We differentiate ([lO]) w.r.t. Xi and for z, j € {1, . . . , d}, that is 
dtdij9 + w • Vdij9 — Adij9 = j^dijw — diu ■ Vdj9 — dju ■ Vdi9 — diju ■ \/9. 
Applying Proposition [2] with r = 1, and 6 = H to Q = dij9 and / = jjdijw - 
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diU ■ ^djO — djU ■ VdiO — diju ■ V6 and summing over all I < i, j < d yields 

{\v*e\)dz 

-H 



< 



M r {\f\)dz 



+ e-^' [ {\Vf\)dz 

J-H 



/ f-H X 1/2 

/r-H \ 1/2 / r-H \ 1/2 

+ iV / (IVti|2)dz / (|V20|2)dz 



\J ~H / \J~H 

/ ,-H \ 1/4 / ,-H \ 3/4 

+ e-^ n ^(iV^xl^) dz\ U jV'e\'/'') dzj 



We invoke (j53p and (|54p successively and estimate with the help of Young's 
inequality 

r {\V^e\)dz < M r mdz + e-^ r {\Vf\)dz 

J~H J~H J~H 



+H+N^ r {\vu\^)dz+e-^^ r 

J-H J-H 



(|Vu|^)(iz.(56) 



Recalling the symmetry properties of the involved quantities at z = 0, we see 
that we can replace the domain of vertical integration in the above estimate by 
[0,i?]. Observe that the (weak) vertical derivative of / at z G {0,-ff} might 
not be defined in the case where dij = d^. In this case, we may interpret 
J {\dzf\) dz in the sense of a BV norm and use the symmetry property to 
estimate 

{\d,f\)dz < 2(|/|,=ol> + 2(|/|.=/^|) + 2 r{\dj\)dz < r{\d,f\)dz, 
H Jo Jo 

where the last inequality follows from a standard Poincare/trace estimate, 
using the fact that d~^f = for z G {0, H}. 



It remains to estimate the terms on the r.h.s. of (|56p . We start with the 
first integral. By the definition of / and using Holder's inequality, we have 



< 1 

^ H 



H / i-H \^/^ / l-H \ 

{\V^w\)dz+{ {\Vu\^)dz] / {\V^e\^)dz 





/ [-H \ 3/4 , H \ 1/4 

+ n {\v^u\'/^)dz\ {m^)dzj 
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Via Young's and Jensen's inequalities and (|52|) -(|54 p . we further estimate 

^ I "-72„„|2\ 



{\f\)dz < {-I ilV^w^dz 







+ r{\Vu\^)dz+ /''(|V2n|4/3)dz + W 



|[29tfcll7t 

< l + HNu+H + {i.) 

< (lnT//)F + W, (57) 

for some 7 > 0. Here and in the following, we denote by (★) terms that can 
be absorbed in the l.h.s. of (l56l). 



For the V/-term, we first compute 

|v/| < + I'^^ll'^^^l + l^^^ll'^^^l + l^^^ll'^^l- 

H 

We estimate using Holder's inequality 

H / i-H \ 1/2 







+ ( 1^ {\v'u\')dzj ^ {\v'e\')dz 

"H \ 1/2 / \ 1/2 

+ ( / {\V\\'^)dz] / {\ve\'^)dz 







We apply Young's and Jensen's inequalities as above and estimate with the 
help of ([53]), dsn) & (021) 



I (|V/|)dz < i^-J^ {\VM')dzj + {\Vu\')dz 

(■H / f-H \ 1/2 

+ y {\V^u\^) dz + ij (|V\|2)dzJ iVwi/2 + (*). 
To control the terms on the right, we finally invoke (|16p and (|29|) : 
j\\Vf\)dz < -L_Nu'/^ + H' + H + Nu+{*) 

< H' + {*). (58) 

Inserting (j57p and ()58p into (j56p and applying ()29p and ([7]) to the remaining 
terms, we have 

/ (|V^0|)(iz < M{ln^H)H + e-^^H^ + H + N^HNu + e~^^H\ 
Jo 
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Choosing M = AlxiH and N = \nH, and using ([S]) and H ^ 1 again leads to 
the desired estimate ([55]) . 

□ 

We now turn to the 

Proof of Theorem[I\ Step 1. Localization of the heat equation. Let r] = 
r]{z) denote a 2if-periodic smooth function with r](z) = Tj{—z) that has 
the cut-off properties = 1 for z G [—1,1], r/ = for z [—2,2], and 
sup(|?7| + [ry'l + |r/"| + \ri"'\) < 1, where ry' = etc. Multiplying ([iO]) by 77 
yields 

dt{rie) + u ■ V {ri9) - /\{r}e) = —riw + r]'w9-2ri'dj-ri"e. (59) 

H 

Step 2. Bound on the second-order derivatives. The bound on the second 
order derivatives relies on Proposition [3l We define Q = rjO and / = -jjr]w + 
r]'w9 — 2r]'dz0 — r]"6 and choose 6 = 2. Eventually after mollifying and / 
and recalling the symmetry of the functions involved, (j50p yields 

2/'^flM2\ <^ / / f2\ I „ / /|V7„.|2\ 



{\V'ir]e)\')dz < / if') dz + sup \ve\ / {\Vu\')dz. (60) 
Jo ^ Jo 

Also notice that the interpolation inequality (j38|) together with (|41|) yields the 
lower bound 

r{\Vm\')dz < r{\V\v0)\')dz. (61) 
Jo Jo 

To bound the r.h.s. of (j60p . we invoke the definition of rj and compute 

ll4T1lfcH>l 

< + (5,0)2 + 1, 

and the r.h.s. is in vertical direction supported in [0, 2]. Hence, using (|4ip and 
the properties of r] again, ([60|) becomes 

[\\V^'9\^)dz < [\w^)dz+ [\{d,ef)dz+ [\\Vuf)dz + l, 
Jo Jo Jo Jo 

and thus, by ([271), (02]) and ([HD, we obtain 

{\V^9\'^)dz < WH, (62) 



/ 

Jo 



where 7 > is an arbitrary constant dependent only on the space dimension, 
and whose value may change from line to line in the remainder of the proof. 
In view of ([6T]) . the same estimate yields 

C{\V9\^)dz < WH. (63) 
^0 
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For the following steps, it is crucial to remark that the estimates in (|6'2p and 
dMD hold in any strip [0, A)"'"^ x (0, h) for /i ~ 1. 

Step 3. Bound on the third-order derivatives. The bound on the third- 
order derivatives relies on Proposition [2] with r = 4/3. Starting point is (j59p . 
Differentiating this localized equation w.r.t. Xi for some i € {1, . . . ,d} yields 

dtdi{7^e) + u-vdi{rie)- Miir^e) 

= —diTjw + —rjdiW + diVj'wO + rj'wdiO 
H H 

- 2dir]'dze - 27]'didze - dir]"6 - 7]"die - rjdiu ■ VO. (64) 

With Q = di{riO) and / defined by the r.h.s. of this equation, and 5 = 2, 
we may apply Proposition [2l Eventually using a mollification argument and 
using the symmetry of C, "W, and /, (|46p yields 

H \ 3/4 

{\v^di{r]e)\'^/^)dz 



H \ 3/4 

H \ 3/4 



< M {^j\\ftri) dz 



{\Vft'^)dz 







H \ 1/2 

mve)f)dz 







4 sl/2. H \ 1/4 



+ N{ 1^ {\Vu\')dzj ^ {{d,{r,e)r)dz 

rH \ 1/4 / rH \ 1/2 



\vuf)dz^ (^^ {\vdi{ve)\^)dz 



Summing over all i, invoking the interpolation inequalities (|38p &: (|39p together 
with (jUJ and using Young's inequality, the above estimate simplifies to 

1 N 3/4 

(|V30|4/3)dz 
/ 

H x3/4 . .H x3/4 



< {\f\'/')dzj {Nf\'^')dzj (65) 

/ M \ 3/4 / r-f^ \ '^/^ 

+ iV3/2M (iV-up)dzJ +^~^^(y (|Vn|^)(izJ +1.(66) 

The velocity terms in are easily controlled. Indeed, we apply and 
(129]) and get 

/ /.4 \ 3/4 . \ 3/4 

iV3/2f / {\Vu\^)dz] +e~''^(/ (|Vn|4)dzj +1 

< In^H, (67) 
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if we choose N = | In i7 and use H ^ 1 in the last inequaUty. We now turn 
to the estimate of the terms in (j65p . By the definition of / and r] and using 
the maximum principle for the temperature, i.e., (j4ip . and H ^ 1, we have 

I/I ^ Inl + lV-ul + lullVei + lVnllV^I + lVei + lV^^I + l, 

and / is supported in [0,2]. Therefore, applying Jensen's and Young's in- 
equality, we estimate 

H \ 3/4 

{\f\'^')dz 



2 \ 1/2 / /.2 \ 1/2 

2\ 



/ .2 X 1/4 / .2 X 1/2 

+ {\ve\^)dzj {\s/^e\^)dzj +1. 



It remains to invoke (|27j) . (j62|) . (j63|) . and if » 1 to deduce 

.H \ 3/4 



We finally investigate the second term in (|65p . Proceeding as in the previous 
computation, but slightly more crudely, we estimate 



|V/| < \u\ + \u\\Ve\ + \u\\V^9\ + \Vu\ + \Vu\\V9\ + \Vu\\V'9\ 
+ |V\| + |V^n||V6l| + |V6l| + IV^^I + \V^e\ + 1. 

We skip the details of how to estimate the norm of V/ step-by-step. 
Via standard inequalities and the global bounds (p9|) and (jSTj) . we can easily 
derive a bound of the form 

rH \ 3/4 

(|V/|4/3)dzj < (69) 



for some q > 0. Hence, choosing M = qlnH and substituting (|67|) - (|69|) into 
(IMD&dMI) yields 

,1 X 3/4 



J {\V^e\^/^)dz] < In'H. 



Step 4- Bounds on the fourth-order derivatives. The derivation of the 
bound on the fourth-order derivatives proceeds very much along the lines of 
that on the third-order derivatives. We differentiate (j64p w.r.t. Xj for some 
j e {!,..., 4: 

^t^^J{r]9) + u■V^^j{r]9)-A^^j{rJ9) = /, 
where / is such that 

I/I ^ Inp + lV-up + lV^-up + lV^I^ + lV^ep + lV^ei^/^ + l, 
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and 

|v/| < \uf + \Vu\^ + \v\\'^ + \v^u\'^ + \ve\'^ + \v^e\'^ + \v^e\^/^ + \v^9\ + i, 

and / is supported in [0,2]. Eventually after mollification, we may apply 
Proposition [2] with r = 1 and 6 = 2 and thus ()46p yields, recollecting the 
symmetries in z = 0: 

H 

{\v^dij{i]e)\)dz 

< M r{\f\)dz 

Jo 

H 







+ e~^^/ {\Vf\)dz 
Jo 



'0 

r-H \ 1/2 

+ ( / {{d^,{ve)f)dz 







1/2 / rH \ 1/2 

{{d,jm)')dz 



1/4 / rH \ 3/4 



+ iv(^^ {\Vu\^)dz^ 

+ e-^ (^j\\Vu\^)dz^ ' [j\\^d,,(r]B)tl'^)dz 

In the case where dzf is not defined as a weak derivative, we argue as in the 
proof of Lemma [3 Summing over all i,j and estimating with the help of 
(j39|) . ([^T]) . (p6]) and Young's inequality, we obtain 



H 

(|V%(7?^)|)dz 



{\^!\)dz 



i-H i-H 

< M / {\f\)dz + e-^' / 
Jo Jo 

+ [ {\Vu\^) dz + e-^^ 
Jo Jo 



4 

4 



ur)dz + l. 

Local bounds on the V^ti-terms can be obtained via the Ehrling estimate 

2„.|2\j^ ^ I /iv7„,|2\ I / /iv73„,|2 



{\V'u\')dz < / {\Vu\')dz+ / {\V^u\')dz. 
Jo Jo 

Proceeding as in the previous steps and observing that the third-order deriva- 
tives of u can be controlled via (|16|) and ([H]), we conclude that 

[\\V^9\)dz < MW H + e-^^H'i + N^\n' H + e~^^H^ + 1. 
Jo 

We choose M = qlnH, N = | In and deduce 



[\\V'e\)dz < In-'H. 
Jo 



□ 
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We finally turn to the proofs of Propositions [2] and [31 The proof of the 
latter one in standard. We display the easy argument for the convenience of 
the reader: 

Proof of Proposition\^ We differentiate ()47p w.r.t. the Xi coordinate, e.g., 

dtdiC + u ■ vdiC - Ad,c = dif - diu ■ vc. 

We multiply this equation by diC and integrate in space and time: 

rH f-H f-H 

{diQdiU ■ VC) dz. 



[ {diCAdiOdz < [ {d,CdJ)dz- [ 

J~H J-H 



'-H 

Here we have used that {dt{di(^)'^) > due to the long-time averages. More- 
over, the transport term drops out because of (|48p . Integration by parts and 
ED yield 



H 

2 



H 



pH r-H 

/ (afc/) dz + / (va,c • diuQ) dz. 

J-H J~H 



It remains to apply the Holder and Young inequalities and sum over 1 < i < d 
to obtain (EOl) via (IMl). □ 



The proof of Proposition [2] requires some preparation. We first provide 
new notation. For a Schwartz function ^l'{x) in W^, we define its Fourier 
transform TTp{q), q G W^, via 



J^^(q) ■= / tl;{x)e"i-'' dx. 

For a periodic function ((x) with cell of periodicity Q = [0, ai] x • • • x [0, a^] 
its Fourier transform T(^(q), q € 27r/aiZ x • • • x 27r/adZ, is defined by 



H{q) ■■= J^J ax)e''-'dx. 



Next, we derive an energy-estimate for narrow-banded solutions to convection- 
diffusion equations. 

Lemma 8. Let u, and f be periodic functions with vertical period 2H 
satisfying 

dtC + u-VC-AC = f. (70) 
Assume that C is narrow-banded in Fourier space in the sense that 

H{q) = for all q B„{qo) (71) 

for some e~^ < [gol < e and < o" ^ 1. For 1 < r < oo it holds 

''{\Cndz< r {\f\^)dz+ r {\V-u\mdz. (72) 

-H J-H J-H 
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Proof of Lemma\^ Let A[s) be a smooth approximation of 

A{s) = -\s\'. 
r 

In view of ([TOD, we have that dtA{C) + u ■ VA{C) - AC^'(C) = fA'{C), and 
thus after integration, integration by parts, and time averaging: 

- {ACA'iO)dz < r {{V-u)A{0)dz+ [" {fA'iO)dz. 

J~H J-H J~H 

Notice that the term involving the time derivative drops out in the time 
average, since {dtA{C,)) > . Now we may carry out the approximation for 
A. Obviously, the assertion in (|72|) follows from 

{ACA'iC))dz = - {AC{signC)\Cr')dz > (KH dz (73) 
and 

{fA'{0)dz = {f{signO\Cr')dz 

/ rH \ lA / i-H \ {r-l)/r 

< / midz] / (icDdz 



-H J \J~H 

via Young's inequality (if r > 1). 

In our argumentation for ()73p , we follow |151 Lemma 1] (or |141 Proposition 
2]). We first show that 

""(l-AC-koPCDd^ < r {\Cndz. (74) 

-H J-H 

For this purpose, we select a Schwartz function ip that satisfies {J-ijj){q) = 1 
for If^l < 1 and define 

An easy calculation shows that {J^'ijja){q) = 1 for g € i?o-(5o)- Therefore, by the 
narrow-bandedness assumption (|7ip . ip^ leaves C invariant under convolution, 
i.e., 

c = V<7 * C- 

It follows that 

-AC-|goPC = (-AVa - |<?oP^<x) * C, 
and thus, applying the convolution estimate yields 

H / /• \r rH 

2^ir\ J. ^ / / I A„;, l„ |2 



{\-AC-\qo\'Cndz<{ \-A^^-\qo\'i^^\dx\ / {\Cndz. 

H KJtR.'' / J-H 

For dZH), we have to show that 

I - Aipa- - \qo\'^'ipa\ dx < a. 
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Indeed, since — A(e ^io-^"j = [^Qpg m-^^ -^q obtain 

so that, because of \qo\ < e and o" < 1 and since ip is a Schwartz function, we 
have 

/ \- Aij„-\qo\'^i!^\dx < a\qo\ {Vipldx + a"^ \V^TP\dx < a. 
It remains to argue how ()74p impUes ()73p . It is 

- {{signOlCr' AC) dz 

= r ((signOlCr-'kol'C) dz + ((signOlCr-^-AC - koPO) dz 

J-H J-H 

> \q,? r midz 

J-H 

H \ {r-l)/r . H \ 1/^ 

^midz^ i^J J-AC-\qo\\ndzj 

> {\q,\^-Ca) {\Cndz 

J-H 

for some constant C > 0. We choose a small enough, such that 2Ca < 
and obtain (1731) because of lonl > □ 



Another tool in the proof of Proposition [2] are the following commutator 
estimates: 

Lemma 9. Let [u-,<j)*]^ denote the commutator of the operations "multipli- 
cation with u" and "convolution with (p", that is, 

[u-, = u ■ * (p) — {u ■ ^) * (p. 

Then we have the estimates 

,-H \ 

< ( f \V^\\i\d±+ [ \cp\dx] 

X i^J^jcr^^P-^'^) dzj i^J jVun dzj , (75) 



and 



{\[u.,(b*]vcndz^ ' 

H 

\<pi^x\ dx 

H \ {p-'^)/{-'-p) / pH \ 1/P 



J ^(|VCr-^/(P-'^)) dzj i^J jVu\n dzj , (76) 

for any 1 < r < oo and 1 < p < oo. 
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Proof of Lemma We only prove the first estimate. The second one is shown 
similarly. We start with a pointwise statement: 

{u ■ (VC * 0) - (n • VC) * (p) (x) 

(/>(x)Vx-C(a; — x) • {u{x) — u{x — x)) dx 

Vx(t>{x) ■ ({x — x){u{x) — u{x — x)) dx 

+ / (l){x)C{x — x){'Vx ■ u){x — x) dx 

< / |V</)(5;)||C(a^ — / \Vu{x — sx)\ dsdx 

Jw' Jo 

+ \(l){x)\\C{x - x)\\{S/ ■ u){x - x)\dx. 
This implies by the triangle inequality 

{\[u.,(b*]s/cndz) ' 

■H J 



|V<^(x)||C(- -x)\\x\ [ |Vn(- -sx)\dsdx] \ dz 

H \ VJRrf Jo ) I . 

+ { l"" (( [ \Hma - -mi"^ -^K - -x)\dxY) dz^ 



< 



-H \ KJ-R.'i 

We consider the first integral only. The estimate of the second one is obtained 
in an analogous way. We have via Holder inequality and Fubini's Theorem 

cH 



H 



r-l 



(^j |V(/)(i)||C(- -x)\\x\ j |Vu(- - sx)\dsd3^ ^ dz 

If \V(t){x)\\x\\(:{- -x)Y f \Vu{- - s£)Y dsdx\ dz 

H \j]R'* Jo I 

)r— 1 

x/ \V^(x)\\x\l [ {\Ci ■ -x)\'\\/ui- - s3:)\'') dzdsdx. 

Ju'i Jo J~H 

We use the Holder inequality again and the translation invariance to estimate 
the inner integrals: 

(|C(- -x)r|Vu(- -sx)Y) dz 

H 



/ pH \ KP-i IIP / nH 

< iyj ^(ic(- -x)r/(f-'^))(izj i^j j^u{ - -sx)\nd, 

H \ {p-r)/p / N 

{\(^\^P/iP-^))dz] U {\\7u\P)dz 
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We combine the above estimates and get 

(^(^j |V(/>(x)||C(- -x)\\x\j^ |Vn(- -sx)\dsdx^ ^ dz 
< U \V^\\x\dx] U {\Wu\P)dzj U {\Cr^^P''-'^)dz] 

□ 

We are now in the position to prove Proposition [21 

Proof of Proposition O The proof reUes on Besov spaces rather than Sobolev 
spaces. We introduce a (non-dyadic) Littlewood-Paley decomposition. Let 
{4>e{x)}i£'z be a family of Schwartz functions with the following properties: 

{T(j)o){q) / only for q with \q\ G (e"\e), (77) 
{T())e){q) = (TMie-^q) ior all£ andq, (78) 
Y,iJ'^i)(^) = lforallg/0. (79) 

We refer to [21 6.1.7 Lemma] for a construction with dyadic blocks, which 
easily adapts to the non-dyadic case. The Littlewood-Paley decomposition 
{<fi}i£'Z of a periodic function (/? is defined via 

ife = (1)1* if. 

Notice that (|77|) - ([79]) imply that Lp = {ip) + ^e- We also introduce the 

low-pass filter 

e<-i 

We now turn to the proof of (|46p . We invoke the triangle inequality to split 
and estimate the integral on the l.h.s. according to 



5 \ 1/^ 



< (|^(|V%<r)dz) (80) 

+ E ^(i^'c^n^^J (81) 

+ E / (iVV)d-- • (82) 



We start with the low frequency part, i.e., (|80l) . It follows from ([77|) - ([79]) 
that 

Co< = * c, 
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where (p^ denotes a Schwartz function with {Tcj)^)(q) = for jg] > e ^. This 
imphes in particular 

J \(j)Q\dx < oo 

and 

l-^(V'(<^o<*C))(g)l = \q\'\H4>^*0{q)\ < TOo<*C)(g)l, 
so that by Plancherel's Theorem and the convolution estimate 

H J-H 
J-H 

< / {e)dz. 

We conclude that 

r5 / rS \ ^'/S 



{\v\<ndz < s('-^-y^ ^1 ^{\v\<\') dz^ 

/ f-H \^/2 

< s('~^y'i^j ^{e)dzj . (83) 



We now address the intermediate frequencies, i.e., (jSip . Our goal is the 
following estimate: 

N-l / i-H \ 



/ [-25 \l/2 / N (2"r)/2r 

+ N{j{\Vu\^)dz\ f y ^(|C|'"/^'~"V^ ) ■ (84) 

Because of (jl5]) . equation ([I3|) is localized in the strip [—5,(5] x R'^"^, so that 
we may replace u in (|43p by u = ryu, where = ri[z) denotes a 2i7-periodic 
extension of a smooth cut-off function satisfying rj = 1 for z € [—6, 6], rj = for 
2; ^ [—25,2(5], sup \r]\ < 1, and sup < ^. We remark that the convection- 
diffusion equation (j43p is invariant under the scaling 

x = e-^x, t = e-^H, C = e"^^C, u = e% f = f, 

i.e., it holds 

diC + ^ • vc - AC = / (85) 

in the rescaled domain of height H = e^H. In a first step, we apply Lemma[8]to 
this rescaled equation. Given some < (7 <C 1 as in the hypothesis of Lemma 
El we select a finite number of open balls Ba{Qj)i<j<j that cover the annulus 
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{e~^ < \q\ < e} and select a family of Schwartz functions {'ipcT,qj}i<j<J that 
form a partition of unity subordinate to that covering. We apply tpa,qj * (po* 
to ([85]) and write the resulting equation as 

^tUqj + ^ • '^C<T,gj - ^Ca,q, = fa,qj + [u-,^a,q, * M^C, (86) 

where subscript ^''CT,qf indicates the convolution with V'o-.gj * <?^o and 
denotes the commutator of the operations "multiplication with u" and "con- 
volution with (/)", that is, 

[u-, (/)*].^ = n • *(/)) — (n • ^) * (/). 
Since Ccr,^^ is narrow-banded in Fourier space in the sense of (j7ip . 

{na,q,){q) = (7-^.,,,)(<?)(-F0o)(g)(J-C)(g) = for ah q B„{qj), 
we infer from (j72p that 



< 



r \ 



H 

Applying the Holder inequality and |V • tt| < |Vn| to the third term 

r(iv-^iic.,,,r)' 

(|Vwp)dz 



< 



dz 



H I 

/ H \ (2-0/(2r) / ^ X (r-l)/r 

the commutator estimate (|75p with p = 2 to the second term 

H 

< 



|V^<x,g^. * (l)o\\x\ dx+ iV'a,?, * M dx 



H / \J-H 



and observing that 

/ W'4^a,qj * 4>o\\x\dx + / iV'o-.gj * f^ol c^a; < oo, (87) 
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since 4'a,qj * <^o is a Schwartz function, the above estimate turns into 



r\ 

.I'^o-.q,- I 

H 

H 



H 



+ I / (|v^P)di 



In view of (j87p . we may use the convolution estimate and Young's inequahty 
(in the case r > 1) to reduce this estimate to 

''A\L,,Mdz 

-H 

< y ^(i/.,,,r)dz+fy ^(jwp)dzj fy^dcr/^^-'^v^j 

By a standard covering argument, we convert the above micro- local inequality 
from the finite number of open balls {Bcj{qj)}i<j<j to a local inequality on 
the annulus {e~^ < \q\ < e}, so that we may replace Ca-,qo by Co and fa,qo by 
/o. Now, we scale back to the original variables. Since (I78p translates into 
(j)i{x) = e"^^ 4>o{e^ x) , it holds Co = e^^Q and fo = fi, and thus 

^2re (^\Q\r)dz 
J-H 

< J Jfendz+i^J ^{\Vu\')dzj jCr/^'~'^)dzj 
Since u = r]u, we have 

H 1-25 1 1-25 p25 

{\Vu\')dz< {\Vu\')dz + -^ {H')dz < {\Vu\')dz. 

■H J-25 J-25 J-25 

The last estimate is due to the Poincare inequality which we may apply thanks 
to the homogeneous boundary conditions of u at z = 0. Now, (j84p immediately 
follows from the well-know fact that 

f" {\V%,\^)dz ~ e^^' {UP)dz, (88) 

for any 1 < p < oo and any s € H, and summing over < ^ < — 1. We 
display the argument for (|88p for the convenience of the reader. As we will 
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not consider applications with fractional or negative derivatives, we restrict to 
the case s G N. Because of ([77|) -([79 |) . we have 1 = X]£'=^-i £ ™ 
support of so that J^cpi = Yl,i'=i-i,t,i+i-^4'i'^4'i- Hence Ylie=i-i,t,i+i't'i-' 
leaves invariant under convolution. We deduce 



so that 

rH rH 



P rH 



i'=i-l,i/~\~l 

< e^P^rHdndz, 



H 



where we have use that 



/ \V''(l)Q\dx < oo. 



For the opposite inequality, observe that because of Q = '^ii^£_i £ £_^_i (pi' * Ci 
we have 



for any 1 < j < d. Thus 
r{\Cindz<( Y I \^^'{^:F<Pi^\d^ r {\dM'')dz. 



We easily compute that 



IT--' 



This proves (|88p for s = 1. The case s > 2 follows by iteration. 

Finally, we consider the high frequency part, i.e., ()82p . We show that 



< Y(rjM")d^) 



/ pH \ 1/4 / X (4-r)/4r 

+ e"^fy (|V'u|4)dzJ U (iVCI^^/^^-^'V-^J • (89) 
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< 



Our treatment is slightly different from the one for the intermediate frequen- 
cies. Starting point is again equation (j86p with u replaced by u: 

diCa,q, + U ■ VCa,qj " Ma,q, = fa,qj + '(/'a.g, * 0O*]VC. 

Using the similar arguments as for the intermediate frequencies but with the 
commutator estimate (j76p and p = 4 instead of (j75p . we arrive at 

" {\v\,ndz 

-H 

/ l-H \ y i-H \ (4-r-)/4 

Notice that the analogous estimate to (f72|) simplifies because of (jM]). Now, 
()89p follows from 



Compared to (|^6|) . it remains to show that 



< Mi^J Jfndzj +e-'^ jVfndzj . (90) 
We split the sum on the l.h.s. of (|9Up according to 

l/r 



Y^lTmidz) 

E (/"(iAr><*^)"^E(/ 



lA / rH \ l/*" 



0<£<A/ 

Thanks to the convolution estimate and 



H 



I \(j)t\dx = / 



6o|da^ < oo, 
we have for the low frequency part 

Q<i<M-l ^•'-^ ^ \J-H 
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and for the high frequencies 

/ f-H \ !/'• (|88} 

< 
< 

We combine the estimates for the high 

ini). 
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